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ABSTRACT 

Wide-field optical imaging surveys will contain tens of thousands of new strong gravi- 
tational lenses. Some of these will have new and unusual image configurations, and so will 
enable new applications: for example, systems with high image multiplicity will allow more 
detailed study of galaxy and group mass distributions, while high magnification is needed to 
super-resolve the faintest objects in the high redshift universe. Inspired by a set of six un- 
usual lens systems [including five selected from the Sloan Lens ACS (SLACS) and Strong 
Lensing Legacy (SL2S) surveys, plus the cluster Abell 1703], we consider several types of 
multi-component, physically-motivated lens potentials, and use the ray-tracing code GLAM- 
ROC to predict exotic image configurations. We also investigate the effects of galaxy source 
profile and size, and use realistic sources to predict observable magnifications and estimate 
very approximate relative cross-sections. We find that lens galaxies with misaligned disks and 
bulges produce swallowtail and butterfly catastrophes, observable as "broken" Einstein rings. 
Binary or merging galaxies show elliptic umbilic catastrophes, leading to an unusual Y-shaped 
configuration of 4 merging images. While not the maximum magnification configuration pos- 
sible, it offers the possibility of mapping the local small-scale mass distribution. We estimate 
the approximate abundance of each of these exotic galaxy-scale lenses to be ^ 1 per all-sky 
survey. In higher mass systems, a wide range of caustic structures are expected, as already 
seen in many cluster lens systems. We interpret the central ring and its counter-image in Abell 
1703 as a "hyperbolic umbilic" configuration, with total magnification ^ 100 (depending on 
source size). The abundance of such configurations is also estimated to be ~ 1 per all-sky 
survey. 
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1 INTRODUCTION 

Strong gravitational lenses, by producing multiple images and high 
magnification, can teach us much about galaxies, how they form 
and how they evolve. The numerous constraints allow us to de- 
termine the mass distribution of galaxies, groups and clusters, in- 
cluding invisible dark substructure. On the other hand, strong lens 
systems, by the high magnification they provide, can act as cosmic 
telescopes, leading to far more observed flux and thereby opening 
unique windows into the early universe. 

Until recen tly, apart from a very few exceptions (e.g. 
iRusin et al.l200lh . the lens configurations observed have been Ein- 
stein ring, double and quadruple image configurations. As the num- 
ber of discovered lenses increases, these configurations will be- 
come commonplace and truly exotic ones should come into view. 
Indeed, several recent strong lens surveys have proved to be just big 
enough to contain the first samples of complex galaxy-scale and 



group-scale strong lenses capable of producing exotic high mag- 
nification image patterns. One such survey is the Sloan Lens ACS 
survey (SLACS, iBolton et alJboodbOOSh , which uses the SDSS 
spectroscopic survey to identify aligrmients of relatively nearby 
massive elliptical galaxies with star-forming galaxies lying behind 
them, and then the Advanced Camera for Surveys (ACS) onboard 
the Hubble Space Telescope (HST) to provide a high resolution im- 
age to confirm that multiple-imaging i s taking place. Anoth er is the 
Strong Lens Leg acy Survey (SL2S, ICabanac et a"i]|2007h , which 
uses the 125-square degree, multi-filter CFHT Legacy Survey im- 
ages to find rings and arcs by their colors and morphology. To- 
gether, these surveys have discovered more than 100 new gravita- 
tional lenses. 

Moreover, strong lensing science is about to enter an excit- 
ing new phase: very wide field imaging surveys are planned for 
the next decade that will contain thousands of new lenses for us to 
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study and expl oit. For example, the Large Synoptic Survey Tele- 
scope (LSST, llvezic et alJIlOOSl) will provide deep 27th AB 
magnitude), high cadence imaging of 20,000 square degrees of the 
Southern sky in 6 optical filters, with median se eing 0.7 arcseconds. 
Similarly, the JDEM and Euclid jLaureiis et al. 2008) space obser- 
vatories are planned to image a similar amount of sky at higher 
resolution (0.1-0.2 arcseconds), extending t he wavelength cover- 
age into the near infra-red (the SNAP design, Alderi ng et al.ll2004l 
with its weak lensing-optimised small pixels and large field of view, 
would be particularly effective at finding strong lenses in very large 
numbers). With such panoramic surveys on the drawing board, we 
may begin to consider objects that occur with number densities on 
the sky of order 10~* per square degree or greater. 

In order to understand the contents of these forthcoming 
surveys, it is useful to have an overview of the configurations, 
magnifications and abundances we can expect from exotic grav- 
itational lenses, namely those with complex mass distributions 
and fortuitously-placed sources. Several attempts have already 
been made in this direction, generally focusing on individual 
cases. This field of research is centred on the notion of criti- 
cal points ("higher-order catastrophes" and "calamities") in the 
3-dimensional space behind a massive object. These points are 
notable for the transition in image multiplicity that takes place 
there, and have been studied in a number of different situa- 
tions by various authors in previous, rather theore t ical works (e.g. 
Blandford & Narayan 1986; Kassiola et al.' 'l992|; | Kormann et al.l 
1994 ; Keeton et al. 2000; Evans & Witt 2001.). 

Here, we aim to build on this research to provide an atlas 
to illustrate the most likely observable exotic lenses. Anticipating 
the rarity of exotic lenses, we focus on the most readily available 
sources, the faint blue galaxies {e.g. lEllislll997h , and investigate 
the effect of their light profile and size on the images produced. 
Our intention is to improve our intuition of the factors that lead to 
exotic image configurations. What are the magnifications we can 
reach? Can we make a first attempt at estimating the abundances of 
lenses showing higher-order catastrophes? Throughout this paper, 
we emphasize the critical points associated with plausible physical 
models, and try to remain in close contact with the observations. 

To this end, we select six complex lenses (three from the 
SLACS survey, two from the SL2S surve y, and the galaxy clus- 
ter Abell 1703 from lLimousin et aill2008bl) to motivate and inspire 
an atlas of exotic gravitational lenses. By making qualitative mod- 
els of these systems we identify points in space where, if a source 
were placed there, an exotic lens configuration would appear. We 
then extend our analysis to explore the space of lens model pa- 
rameters: for those systems showing interesting caustic structure, 
that is to say presenting higher-order catastrophes, we study the 
magnifications we can reach, and make very rough estimates of the 
cross-sections and abundances of such lenses. 

This paper is organised as follows. In Section|2]we provide a 
review of the basic theory of gravitational lenses and their singulari- 
ties, and then describe in Section[3]our methodology for estimating 
the source-dependent magnifications and cross-sections of gravi- 
tational lenses. We then present, in Section |4l our sample of un- 
usual lenses, explaining why we find them interesting. Suitably in- 
spired, we then move on the atlas itself, focusing first on simple, i.e. 
galaxy-scale, lens models (Section[5]l, and then on more complex, 
i.e. group-scale and cluster-scale lenses (Section[6]l. Finally, we use 
the cross-sections calculated in Sections [5] and [6] to make the first 
(albeit crude) estimates of the abundance of exotic lenses for the 
particular cases inspired by our targets in Section|7] We present our 
conclusions in Section[8] When calculating distances we assume a 
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Figur e 1. Schema of a typ ical, single plane, gravitational lens system - 
from Schneider et alj j2006l Figure 12, p. 20, reproduced with permission 
©ESO). 



general-relativistic Friedmann-Robertson- Walker (FRW) cosmol- 
ogy with matter-density parameter Q,m ~ 0.3, vacuum energy- 
density parameter Qa ~ 0.7, and Hubble parameter Ho ~ 70 
km s"^ Mpc"^. 



2 GRAVITATIONAL LENSING THEORY 

This section is divided into two parts. We first briefly review the 
relevant basic equation s in lensing theory, following the notation of 
ISchneidere7ai] ( l2006l) . We then introduce th e qualitative elements 
of lens singularity theory l lPetters et alj200l1) relevant to this study. 



2.1 Basics 



We reproduce in Fig.[T]a sketch, from lSchneider et al] ( |2006[) . of a 
typical gravitational lens system. In short, a lens system is a mass 
concentration at redshift Zd deflecting the light rays from a source 
at redshift z^. The source and lens planes are defined as planes per- 
pendicular to the optical axis. 

The lens equation relates the position of the source to its ob- 
served position on the sky; in the notation of the sketch, /3 — 
— a{0) where (3 = (/3i,/32) denotes the angular position of 
the source, 9 = (61,62) the angular position of the image and 
a{ff) — ^jy^a{Dd0) is the scaled deflection angle. D^, Da, and 
Dds are the angular diameter distances to the lens, to the source, 
and between the lens and the source respectively. 

The dimensionless surface mass density or convergence is de- 
fined as 



E(I>d6») 



with Ecrit 



Da 



(1) 



47rG DdDas 

where the critical surface mass density Ecrit delimits the region 
producing multiple images (within which (k) ^ 1). The critical 
surface mass density is, therefore, a characteristic value for the sur- 
face mass densi ty which is the dividin g line between "weak" and 
"strong" lenses dSchneider et ai]|200d p. 21). 

Gravitational lensing conserves surface brightness, but the 



Exotic Gravitational Lenses 3 



shape and size of the images differs from those of the source. The 
distortion of the images is described by the Jacobian matrix, 



^^^^ - de - [^'^ ~ W- ) 



1 — K - 7i —72 
—72 1 — K + 71 



(2) 
(3) 



where we have introduced the components of the shear 7 = 71 + 
172. Both the shear and the convergence can be written as combi- 
nations of partial second derivatives of a tensing potential if): 



0^.11 + i',22) 



71 = - -^,22), 72 = '/',i2- 



(4) 
(5) 



The inverse of the Jacobian, M(0) = J ^, is called the magnifica- 
tion matrix, and the magnification at a point 60 within the image is 
|^t(0o)|, where ^ = det M. 

The total magnification of a point source at position 6s/3 is 
given by the sum of the magnifications over all its images, 



(6) 



The magnification of real sources with finite extent is given by the 
weighted mean of /ip over the source area. 



/d2/3/(^)(/3)Mp(/3) 



[/d2/3/W(/3)] ' 
where (/3) is the surface brightness profile of the source. 



(7) 



2.2 Singularity tlieory and caustic metamorphoses in 
gravitational lensing 

The critical curves are the smooth locii in the image plane on which 
the Jacobian vanishes and the point magnification is formally in- 
finite. The caustics are the corresponding curves, not necessarily 
smooth, obtained by mapping the critical curves into the source 
plane via the lens equation. 

A typical caustic presents cusp points connected by fold lines, 
which are the g eneric singularities associated with strong lenses 
jKassiola et ( Il992h ). By this we mean that these two are sta- 
ble in the source plane, i.e. the fold and cusp are present in strong 
lenses for all Zb > za_. By considering a continuous range of source- 
planes, the folds can be thought of as surfaces in three-dimensional 
space extending behind the lens, while cusps are ridgelines on these 
surfaces. Indeed, this is true of all caustics: they are best thought 
of as three-dimensional structures (surfaces) lying behind the lens, 
which are sliced at a given source redshift and hence renderable 
as lines on a 2-dimensional plot. We will show representations of 
such three-dimensional caustics later. Other singularities can ex- 
ist but they are not stable, in the sense that they form at a specific 
source redshift, or in a narrow range of source redshifts: they repre- 
sent single points in three-dimensional space. The unstable nature 
of these point-like singularities can be used to put str ong constrain ts 
on lens models and can lead to high magnifications teaglallOOlh - 
it is these exotic lenses that are the subject of this paper. 

Excellent introductions to gravitation al lenses and their crit- 
ical points can be found in the books by ISchneider et alj j"l992l , 
where chapter 6 is particularly relevant) and Fetters et alj ( 12001 ): 
here we provide a brief overview. We will consider at various points 



the following critical points (known as "calamities" or "catastro- 
phes"): lips, beak-to-beak, swallowtail, elliptic umbilic and hy- 
perbolic umbilic. With these singularities, we associate metamor- 
phoses: the beak-to-beak and lips calamities, and the swallowtail 
catastrophe, mark the transition (as the source redshift varies) from 
zero to two cusps, the elliptic umbilic catastrophe from three to 
three cusps, and the hyperbolic umbilic catastro phe involves the ex- 
change of a cusp between two caustic lines dSchneider et alJ[T992l 
section 6.3). These are the five types of caustic metamorphosis that 
arise generically from taking planar slices of a caustic surface in 
three dimensions. As well as these, we will study the butterfly meta- 
morphosis, which involves a transition from one to three cusps. 
It can be constructed from combinations of swallowtail metamor- 
phoses; lenses that show butterfly metamorphoses will have high 
image multiplicities, hence our consideration of them as "exotic" 
systems. IPetters et al. U2OOII section 9.5) give an extensive discus- 
sion of these caustic metamorphoses, and we encourage the reader 
to take advantage of this. 

In Fig, p and F ig. [3] which are reproduced and extended from 
IPetters et al.l 1 200 ih . we sketch these critical points and their re- 
lated metamorphoses, showing the changes in image multiplicity 
at the transition. The features in the caustics after the swallowtail 
and butterfly metamorphoses are also referred to as, respectively, 
swallowtails and butterflies after their respective shapes. A caus- 
tic shrinking to an elliptic umbilic catastrophe and then appe aring 
again is called a "deltoid" caustic jSuvu & Blandfordll2006l) . We 
note that varying source redshift is not the only way to bring about 
a metamorphosis: the independent variable could also be the lens 
component separation in a binary lens, the isopotential ellipticity, 
the density profile slope and so on. In this paper we aim to illustrate 
these metamorphoses in plausible physical situations. 

A number of theoretical studies of these c r itical points have 
already been performed. .Blandford & NaravanI ( 11986 ) undertook 
the first major study of complex gravitational lens caustics, and 
provided a classifica tion of gra vitational lens images. Following 
this initial foray, K assiola et alj (1992) studied the lips and beak- 
to-beak "calam i ties," applied to the long straight arc in Abell 2390. 
iKormann et presented an analytical treatment of several 

isothermal elliptical gravitational lens models, including the singu- 
lar isothermal ellipsoid (SIE) and the non-singular isothermal el- 
lipsoid (NIE). Keeton et al. (2000) focused on the most common 
strong lenses, namely isother mal elliptical density galaxies in the 
presence of tidal perturbations. Ie vans & Wit3 ( l200lh presented, via 
the use of boxiness and skewness in the lens isopotentials, the for- 
mation of_swallowtail and butterfly catastrophes, while most re- 
cently IShin & EvansI 12008b) treated the case of binary galaxies 
case, outlining in particular the formation of the elliptic umbilic 
catastrophe. 

The mathematica l definitions of th e various catastrophes 
(given in chapter 9 of IPetters et alj|200ll and the above papers, 
along with some criteria for these to occur), involve various spa- 
tial derivatives of the Jacobian, as well as the magnification itself. 
This is the significance of the "higher-order" nature of the catastro- 
phes: they are points where the lens potential and its derivatives are 
very specifically constrained. This fact provides strong motivation 
for their study - lenses showing image configurations characteris- 
tic of higher-order catastrophes should allow detailed mapping of 
the lens local mass distribution to the images. In addition, the large 
changes in image multiplicity indicate that these catastrophe points 
are potentially associated with very high image magnifications, per- 
haps making such lenses very useful as cosmic telescopes. 
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Figure 2. Caustic metamorphoses and their lensed image multiplicities. The 
points determine the critical points. The numbers N, N + 2, N + A indicate 
how many lensed images are produced of a hght source lying in the given 
region. The caustic curves are oriented so such that, locally, light sources 
to the left genetically have t wo more lensed ima ges than those to the right. 
Adapted from Figure 9. 19 in IPettersetalJtoOlL p. 382). With kind permis- 
sion of Springer Science+Business Media. 
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Figure 3. A butterfly caustic metamorphosis, i ncluding a swallowt ail as an 
intermediate step. Adapted from Figure 9.17 in lPetters et all i200lL p. 379). 
With kind permission of Springer Science+Business Media. 



3 SOURCE SIZES AND OBSERVABLE MAGNIFICATION 

We are interested in the most likely instances of exotic lensing 
by systems close to catastrophe points. The most common strong 
lenses are galaxies lensed by galaxies, and so it is on these systems 
that we focus our attention. The observability of strong lensing de- 
pends quite strongly on the magnification induced - and this in turn 
depends on the source size, as indicated in Equation[71 In this sec- 
tion we investigate this effect. 




Figure 4. Upper panel: critical curves (left) and caustics (right) of an 
elliptically-symmetric galaxy modelled by an NIE profile (appendix). The 
source can be seen lying close to a cusp, giving three merging images, 
a fainter counter-image and a de-magnified central image. Lower panel: 
source-plane total magnification maps convolved with four different source 
types: point source, de Vaucouleurs bulge, exponential disk, Gaussian. 
Contours of equal magnification are shown on all three convolved maps, 
overlaid on the arbitrary greyscale images: note the higher magnifications 
reached with the peakiest de Vaucouleurs profile. 



We compute maps of the total magnification in the source 
plane, /ip(/3) in Equation[6l using the ray-tracing code GLAMROC 
(described in the appendix). In order to study the magnification of 
realistic extended sources, we convolve these magnification maps 
by images of plausible sources, described by elliptically symmetric 
surface brightness distributions with Sersic - Gaussian, exponential 
or de Vaucouleurs - profiles (see the appendix for more details). We 
then use these to quantify how likely a given magnification is to ap- 
pear by calculating: a) the fraction of the multiply-imaged source 
plane area that has magnification above some threshold, and con- 
versely b) the magnification at a given fractional cross-sectional 
area. 

We first use these tools to analyze the likely effect of source 
size and profile on the observable magnification. We consider a 
standard lens model, with a non-singular isothermal profile (see 
appendix), and a fiducial source at redshift Zs = 2 with half-light 
radius 2.9kpc (Ferguson et al .l2004h . The lens model and four mag- 
nification maps are presented in Fig.|4] In Fig. |5]we then plot the 
corresponding fractional cross-section area as a function of thresh- 
old total magnification. We first remark that the area of multiple 
imaging corresponds approximately to the area with magnification 
greater than 3, a value we will take as a point of reference in later 
sections. The behaviour of the three source profiles is fairly similar, 
with the peakier profiles giving rise to slightly higher probabilities 
of achieving high magnifications. The two profiles we expect to be 
more representative of faint blue galaxies at high redshift, the expo- 
nential and Gaussian, differ in their fractional cross-sectional areas 
by less than 10%. 

We now turn to influence of the half-light radius on the ob- 
servable magnification, considering only the exponential profile 
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Figure 5. The fractional cross-sectional area (relative to the multiple imag- 
ing area) as a function of threshold total magnification, for a model NIE 
galaxy. The source half-light radius was fixed at 2.9kpc. 





Figure 6. Upper panel: for an exponential profile source, we plot the frac- 
tional cross-sectional area for magnifications greater than various magnifi- 
cation thresholds, as a function of the half-light radius (in kpc, Zg =). 
Lower panel: the corresponding magnification for a given fractional cross- 
sectional area threshold, also as a function of re. 



and keeping the fiducial source redshift of Zs = 2. We plot two 
graphs in Fig.|6] both as a function of half-light radius r^, we show 
the fractional cross-sectional area greater than some magnification 
threshold, and also the magnification corresponding to a given frac- 
tional area threshold. We study an inclusive ra nge of half light-radii, 
from 0.5 kpc to 5 kpc i lMarshall et aP feOOTh give an example of a 
compact source of half-light radius ~0.6 kpc). 

We see that the magnification is, indeed, highly dependent on 



the source size: the smaller the source, the more likely a large 
total magnification is (lower panel). At small magnification, the 
fractional cross-sectional area reaches a maximum: larger sources 
smear out the magnification map more, while smaller sources have 
already resolved the magnification map at that threshold and cannot 
increase the interesting area. 

Clearly the size of the source is more important than the form 
of the source profile. Henceforth we use only one profile, the ex- 
ponential, and study three different sizes: 0.5 kpc, 1.5 kpc and 2.5 
kpc, keeping these constant in redshift for simplicity in the analy- 
sis. At low r edshift, these sizes are somewhat smaller than the ones 
discussed in lFerguson et al.l ( |2004) (and the appendix). One might 
argue that such an arbitrary choice is not unreasonable, given the 
uncertainty and scatter in galaxy sizes at high redshift. The mag- 
nification bias noted here will play some role in determining the 
abundance of lensing events we observe; however, as we discuss 
later, we leave this aspect to further work. In this paper we simply 
illustrate the effects of source size on the observability of various 
lensing effects, noting that 0.5-2.5 kpc represents a range of plau- 
sible source sizes at zs>, 1. 



4 A SAMPLE OF UNUSUAL LENSES 

The majority of the known lenses show double, quad and Einstein 
ring image configurations well reproduced by simple SIE mod- 
els. There are several cases of more complex image configurations 
due to the presence of several sources; there are also a few more 
complex, multi-component lens systems, leading to more complex 
caustic structures. One example of such a complex and interest- 
ing lens is the B1359+15 4 system, compose d of 3 lens galax- 
ies giving rise to 6 images iRusin et alj j200lh . Another example 
is MG2016+112, whose quasar images have flux and astrometric 
anomalies attributable to a small satellite companion (?, and refer- 
ences therein). 

We here show here six unusual lenses that have either been 
recently discovered or studied: three from the SLACS survey 
teolton et alj |2008|) . two from the SL2S surve y ('Caban ac etaU 
l2007h . and the rich cluster Abell 1703 l lLimousin et al. 2008lJ). This 
set presents a wide range of complexity of lens structure. These 
lenses are presented in this section but are not subject of further 
studies in the following sections: their function is to illustrate the 
kinds of complex mass structure that can arise in strong lens sam- 
ples given a large enough imaging survey. All six have extended 
galaxy sources, as expected for the majority of strong lenses de- 
tected in optical image survey data. 

The three SLACS targets are massive galaxies at low redshift 
(zd ~ 0.2), that are not as regular as the majority of the SLACS 
sample jBolton et al. 2008), instead showing additional mass com- 
ponents such as a disk as well as a bulge, or a nearby satellite 
galaxy. The two SL2S targets are more complex still, being com- 
pact groups of galaxies. The presence of several lensing objects 
makes them rather interesting to study and a good starting points for 
our atlas. The cluster Abell 17 03 was first iden t ified a s a strong lens 
in the SDSS imaging data by iHennawi et al 1 ll2008h . and was re- 
cently studied with//5rbv iLimousin et aLlfeOOSbh TAs we will see, 
elliptical clusters such as this can give rise to a particular higher- 
order catastrophe, the hyperbolic umbilic. 
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Figure 7. Three unusual SLACS lenses. HST/ACS images in the F814W filter, approximately 6 arcsec on a side. Left: SDSSJl 100+5329; centre: 
SDSSJ0808+4706; right: SDSSJ0956+5100. In this final image, we denote: 1) the main lens galaxy; 2) the satellite lens galaxy; 3) the counter-image to 
the extended arc system. 



4.1 The three SLACS lenses 

Three unusual SLACS lenses are presented in Fig. [7] 
SDSSJllOO+5329 presents a large circular arc relatively far 
from the lens (~ 1.5"), broken at its centre. We were able to 
model this system with two SIE profiles of different ellipticities 
but with the same orientations, representing a mass distribution 
comprising both a bulge and a disk. While a detailed modeling will 
be presented elsewhere, we take a cue from SDSSJl 100+5329 
and study in Section [5] the catastrophe points present behind 
disk-plus-bulge lenses. 

SDSSJ0808+4706 is a massive elliptical galaxy with a large 
satellite galaxy lying just outside the Einstein ring and distorting 
its shape. SDSSJ0956+5100 also has a satellite galaxy, this time 
inside the Einstein radius (the arcs due to the lensing effect are at 
larger radius than the satellite). The image component labeled 3 in 
Fig. |7] is considered to be due to the lensing effect: the system of 
images to which it belongs is clearly distorted by the satellite. In- 
spired by these two systems, we study below the cata strophes lying 
behind binary lenses, a subject also investigated bv lShin & Evans! 
( l2008ah . 



4.2 The two SL2S lenses 

We c onsider two lenses from the SL2S survey jCabanac et al.l 
l2007h : SL2SJ0859-0345 and SL2SJ1405+5502. In 
SL2SJ1405+5502 we have an interesting binary system leading 
to an asymmetric "Einstein Cross" quad configuration (Figure [8]l. 
Binary systems, as we will see in Secti on [6] below can develo p 
exotic caustic structures (as discussed by IShin & Evand b008al) ). 
We are able to reproduce the image configuration fairly well with 
a simple binary lens model, and present its critical curves and 
caustics in Section|6]below. 

SL2SJ0859— 0345 is more complex, consis ting of 4 or 5 
lens co mponents in a compact group of galaxies ( Limo usin et al.l 
l2008ah . The image configuration is very interesting: an oval-shaped 
system of arcs, with 6 surface brightness peaks visible in the 
high resolution (low signal-to-noise) HSTfWFPCl (Wide-Field and 
Planetary Camera 2) image. It is in this single-filter image that 
the fifth and faintest putative lens galaxy is visible; it is not clear 



whether this is to be considered as a lens galaxy, or a faint central 
image. In Fig.[9]we show the WFPC2 and CFHTLS images, and in 
the central panel a version of the latter with the four brightest lens 
components modelled with GALFIT (in each filter independently) 
and subtracted. The central faint component remains undetected 
following this process, due to its position between the sidelobes of 
the fit residuals. We will improve this modelling elsewhere; here, 
in Section[6] we simply investigate a qualitatively successful 4-lens 
galaxy model of the system, focusing again on the exotic critical 
curve and caustic structures and the critical points therein. 

4.3 Abell 1703 

The strong lensing galax y cluster Abell 1703 (Zd = 0.28, shown in 
Fig. [To] reproduced from'Limo usin et ^1 12008^ exhibits an inter- 
esting feature: a central ring with a counter-image. A naive interpre- 
tation would link this central-quad feature to the two small galaxies 
lying inside it; however, the large image separation across the quad 
and the presence of the spectroscopically-confirmed counter-image 
argue against this. We instead interpret the image configuration as 
being due to a source lying close to a hyperbolic umbilic catas- 
trophe, a critical point anticipated for elliptically-extended clusters 
such as this one ( lKassiolaetalJl992h . 



5 SIMPLE LENS MODELS 

Motivated by our SLACS targets, we study in this section two dif- 
ferent "simple" lens models: a single galaxy composed of both a 
bulge and a disk, and a main galaxy with a companion satellite 
galaxy. In this respect the adjective "simple" could also be read as 
"galaxy-scale." We will investigate some "complex," group-scale 
lens models in the next section. 

Aa single isothermal elliptical galaxy yields at most four visi- 
ble images; higher multiplicities can be obtained when considering 
deviations fro m such profiles, lead ing to more than four images, as 
shown by e.g. jKeeton et all2000l) . Likewise, jEvans & Witll200lh 
showed how boxy or disky lens isodensity contours can give rise to 
six to seven image systems. Here we construct models from sums 
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Figure 9. The complex lens system SL2SJ0859— 0345. Left: g-r composite image from CFHTLS; centre: after subtraction of the four bright lens components; 
right: at higher resolution, with HSTIWVPC2 in the F606W filter. Images are approximately 15 arcsec on a side. 





Figure 8. SJ 1405+5502. found in the CFHTLS data {top, g-r-i composite) 
and followed-up with HSTfWPVCl in Snapshot mode with the F606W filter 
{bottom). Both images are approximately 6 arcsec on a side 



Figure 10. Centeal pa rt of the cluster Abell 17 03, imaged with //ST/ ACS 
and reproduced from iLimousin et al ] j2008bl) . We are interested in the 
bright central ring (images 1.1-1.4), and its counter-image (1.5). 



of mass components to achieve the same effect, associated with the 
swallowtail and butterfly metamorphoses. 



5.1 Galaxies with both bulge and disk components 

5.1.1 Model 

We model a galaxy with both a disc and a bulge using two concen- 
tric NIE models with small core radii, different masses and different 
ellipticities. These two NIE components are to be understood as a 
bulge-plus-halo and a disk-plus-halo components, as sugg ested by 
the "c onspiratorial" results of the SLACS survey jKoopmans et al.l 
l200d) . For brevity we refer to them as simply the bulge and disk 
components. 

We fix the total mass of the galaxy to be that of a fidu- 
cial massive early-type galaxy, with overall velocity dispersion of 
a = 250km s"\ and impose the lens Einstein radius to correspond 
to this value of a. We then divide the mass as Ubuige = 200 and 
f7"disk = 150km s~^ to give a bulge-dominated galaxy with promi- 
nent disk, such as that seen in SDSSJl 100+5329. Integrating the 
convergence, we calculate the ratio of the mass of the bulge compo- 
nent to the total mass, and find / = = 0.64, where Me 
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indicates the mass enclosed within the Einstein radius . This ratio is 
somew hat low for a lens galaxy, but not extremely so: lMoller et al.l 
studied a theoretical population of disky lenses, and found 
that around one quarter of their sample had bulge mass fraction less 
than this. 

The parameters of our double NIE model are summarised as 
follows: 



Type 


a 


e 


Tc (t> 


Bulge 


200 


0.543 


0.1 0.0 


Disk 


150 


0.8 


0.1 4> 



Here the variable 4> represents the relative orientation of the two 
components' major axes, and a is measured in km s~^, in kpc, 
9 in arcseconds and (/> in radians. The ellipticity e refers to the 
mass distribution (see the appendix for more details). Predicted lens 
system appearances, and critical and caustic curves, are shown in 
Fig- EH for three values of the misaligimient angle (f>. 

5.1.2 Critical curves, caustics and image configurations 

The mapping of the outer critical curve to the source plane (right 
panels of Fig. lilt gives quite complex caustic curves, depending 
of the orientation. For = 0, the profile is quasi-elliptical, and the 
inner caustic is a familiar astroid shape. As increases, the astroid 
caustic becomes skewed. At ~ 1, two of the folds "break" in 
to two swallowtail catastrophes (at opposing points). This critical 
value of (f) varies with disk-to-bulge ratio, and the ellipticities of 
the two components: larger ellipticities or disk-to-bulge ratios give 
more pronounced asymmetry, and a smaller is required to break 
the folds. For (j) > (jicrit, we observe two swallowtails present- 
ing two additional cusps each. As <f) approaches n/2, each swal- 
lowtail migrates to a cusp, and produces a butterfly (we can also 
obtain two butterflies with (j) < n/2, if we increase the asym- 
metry). If the asymmetry is sufficiently large, the two butterflies 
overlap, producing a nine-imaging region (eight visible images and 
a de-magnified central image). In the same way, two swallowtails 
can overlap, producing also a seven-imaging region (six images in 
practice). Evans & W itt (2001) and Keeton et al. ( 2000i) give ana- 
lytical solutions to produce such regions by using, respectively, de- 
viations of isophotes (boxiness and skew) and by adding external 
shear corresponding to external perturbations. All we show here is 
that the combination of realistic disk and bulge mass components 
is a physical way to obtain such distorted isodensity contours; the 
disk-to-bulge ratio and their misalignment can be straightforwardly 
estimated using standard galaxy morphology tools. 

Observationally speaking, it is important to note that the var- 
ious critical areas (the swallowtail and butterfly, and the regions 
of overlap between them) are very small, and close to the optical 
axis of the lens. A small source lying in one of these regions will 
appear either as a very large arc (four images merging in the case 
of a swallowtail, five images merging in the case of a butterfly), 
or as a "broken" Einstein ring (in the nine-imaging region cases). 
However, in practice, with realistic source sizes, these image con- 
figurations will appear as Einstein rings, around which the surface 
brightness varies: this is illustrated in Fig.[T2] where we predict op- 
tical images for the same lens models as in Fig.ll llbut with a source 
at Za = 1.2 with — 2 kpc. These fluctuations can tell us some- 
thing about the detailed structure of the large-scale po tential, and 
should be borne in mind when m odelling such rings ( iKoopmani 
l2005l : IVegetti & Koopmanll2009h . 
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Figure 13. Fractional cross-section as a function of total magnification 
produced for our model disk-plus-bulge galaxy, with relative orientation 
= 1.55. 



5.1.3 Magnification 

In order to study the magnification of this model system, we first 
plot in Fig. I13the fractional cross-sectional area as a function of 
the magnification for the almost maximally-misaligned case, — 
1.55. We show three different source sizes (as in Section|3]l. For a 
typical extended source (ro = 2.5 kpc), the total magnification only 
reaches a maximum of around p = 5, no larger than we might ex- 
pect for a simple NIE lens. At fixed source size, Einstein rings have 
similar magnification regardless of whether they have catastrophe- 
induced fluctuations. We can see from Fig. [TT] that when we in- 
crease the misalignment (f>, the astroid caustic area decreases. As 
the source is typically more extended than this astroid caustic, the 
regions of high magnification have a small contribution to the frac- 
tional cross-section. When we decrease the size of the source, the 
magnification shifts to higher values (10-15), but still comparable 
to those obtainable with a simple NIE model. 

We then compute the fractional cross-section of the lens hav- 
ing a magnification higher or equal than a given threshold. Fig. 1141 
as explained in Section [3] With a source size of re = 0.5 kpc, we 
see that for a low magnification threshold (n > 8), the maximum 
is at — 1.25 and not at (f} = 1.57 as we expected. Increasing 
the threshold magnification, the maximum does shift to <jf> = 1.57. 
Increasing the misalignment 0, regions of higher magnification but 
smaller fractional area appear, while at the same time the astroid 
caustic cross-section decreases. By considering higher magnifica- 
tion thresholds, we focus on more the higher magnification regions. 
For example, at tf) — 1.25, the caustic presents two relatively large 
swallowtails, which are regions of medium-high magnification. For 
higher magnification thresholds (e.g. /i > 15), the swallowtail re- 
gions (and indeed the four cusp regions) no longer provide sig- 
nificant cross-section, and the maximum shifts to the butterfly re- 
gions (and the region where they overlap). However, this region 
is correspondingly smaller, and the cross-section decreases to just 
~ 5 X 10~^ of the multiple imaging area. 

These considerations can be used to try and estimate the frac- 
tion of butterfly configurations we can observe. The problem for 
the swallowtails is to find a magnification threshold which isolates 
these regions (and separates them from the effects of the 4 cusps). 
The butterfly regions are more cleanly identified in this way: only 
the place where the two butterflies overlap can give rise to such a 
small fractional area yet high magnification. When estimating but- 
terfly abundance in Sectionl?] we will make use of the approximate 
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Figure 11. Gravitational lensing due to a model bulge-dominated, disky galaxy. Three different bulge-disk misalignment angles are shown: from top to 
bottom, </) = 1-0, 1.3, 1.57 radians. Left to right: a predicted high-resolution optical image; the critical curves and image positions (note the small inner 
critical curve arising from the presence of the small core radius in the NIE profile); the source position relative to the caustics. For the anti-ahgned case, 
<p = 1.57, two butterflies overlap in the centre of the source-plane giving rise to an Einstein ring image system composed of 8 merging images, visible here 
for the case of an almost point-like source. 



fractional cross-sectional area of ~ 5 x 10^'^. For the swallowtails, 
we note that these form over a wider range of misalignment angles, 
perhaps ±0.5 rad from anti-alignment compared to ±0.05 rad for 
butterfly formation: we predict very roughly that we would find ten 
times more swallowtails than butterflies just from their appearing 
for a ten times wider range of disk/bulge misalignments. 



5.2 Galaxies with a satellite 



As illustrated in sectiongjwith the system SDSSJ0808±4706 and 
SDSSJ0956±5100, main galaxies can have small satellites, rela- 
tively close. We expect these satellite to perturb the caustic struc- 
tures, and perhaps give rise to higher order catastrophes (e.g. ?). 
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Figure 12. Gravitational lensing due to a model bulge-dominated, disky galaxy, and a realistic faint blue galaxy source. The same three bulge-disk misalignment 
angles <b as in Fig. [ID are shown: </> = 1.0, 1.3, 1.57 radians, from left to right. The source is at Zs = 1.2 and has a half-light radius of 2 kpc. 




A ^ (radj A ^ {rod) A ^ {rod) 



Figure 14. Fractional cross-section area giving total magnification greater than a given threshold, as a function of disk/bulge misalignment angle. We consider 
our 3 fiducial sources, from left to right: rc = 0.5, 1.5, 2.5 kpc 



5.2.1 Model 

We model a system comparable to SDSSJ0808+4706, i.e. a main 
galaxy with a small satellite, by two NIE profiles (with very small 
cores). The model parameters were set as follows: 



Type 


a 






01 


62 





NIE 


236. 


0.05 


0.05 








1.3 


NIE 


80. 


0.15 


0.01 


-0.735 


1.81 


-0.3 



where a is in km s , rc in kpc, 6 in arcseconds and in radians. 

The redshifts of the lens and source planes were set at Zd = 
0.2195 and = 1.0215 (to match those of SDSSJ0808+4706). 
We arbitrarily chose the velocity dispersion of the satellite to be 
one third that of the main galaxy, and used a double-peaked source 
to better reproduce the arcs seen in SDSSJ0808+4706. 

5.2.2 Critical curves, caustics and image configurations 

The critical curves and caustics for this model system are shown in 
Fig. [15] The two astroid caustics merge to give an astroid caustic 
with six cusps instead of four, a beak-to-beak metamorphosis (see 
section [Z2] l. 

In Fig. [T6l we present the beak-to-beak calamity image con- 
figuration by slightly adjusting the velocity dispersion parameter of 
the satellite from cr = 80 km s^^ to cr = 67.7 km s~^. The image 
configuration of the beak-to-beak configurati on is the merging of 
three images in a straight arc (we refer also to lKassiola et al J 19921 , 
who discuss the beak-to-beak in order to explain the straight arc 
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Figure 15. Critical curves and caustics for our model of SDSSJ0808-f 4706, 
with predicted images overlaid. 



in Abell 2390). The beak-to-beak calamity does not lead to higher 
magnification (compared to the cusp catastrophes), as it is still the 
merging of just three images. 

As the lens potential of the satellite increases, the caustics 
become more complex; such lenses are better classified as "bi- 
nary" and are studied below in Section |6] On the other hand, 
small (and perhaps dark) satellite galaxies are expected to mod- 
ify the lensing flux ratio {e.g. Shin & Evans 2008b): such substruc- 
ture can easily creat e very local catastrophes such as swallowtails 
teradac et al.ll2004ah . which may be observable in high resolution 
ring images, as discussed in the previous section l lKoopman3l2005l : 
IVegetd & KoopmansllIoO^ . 
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Figure 16. Critical curves and caustics for a model main galaxy plus satel- 
lite lens, presenting a beak-to-beal( calamity, and a resulting image config- 
uration. 



6 COMPLEX LENSES 

We expect group and cluster-scale lenses to produce higher 
multiplicities, higher magnifications, and different catastro- 
phes. Motivated by two SL2S targets, SL2SJ1405-I-5502 and 
SL2SJ0859-0345, and the cluster Abell 1703, we first study a bi- 
nary system, and the evolution of its caustics with lens component 
separation and redshift, and the apparition of elliptic umbilic catas- 
trophes. Then we illustrated a four-lens component system based 
on the exotic SL2SJ0859— 0345 lens. Finally, we study the clus- 
ter Abell 1703, which presents an image system characteristic of a 
hyperbolic umbilic catastrophe. 



6.1 Binary lenses 

6.1.1 Model 

Given the asymmetry of the pair of lens galaxies visible in 
SL2SJ1405+5502 (Section [4.2t . we expect the critical curves and 
caustics to be quite complex. We first make a qualitative model of 
the lens, using the following mass components: 



type 


a 


/ 




01 


02 




NIE 
NIE 


250 
260 


0.06 
0.04 


0.1 
0.1 


0.0774 
-0.0774 


-0.55 
0.55 


1.5 
-0.6 



where a is in km s^^, rc in kpc, 9 in arcseconds and (j> in radians. 
The redshifts chosen are = 0.6 and = 1.2, reasonable esti- 
mates for the SL2S objects. As Fig. [TT] shows, the model-predicted 
image configuration matches the observed images fairly well. We 
note the presence of a predicted faint central image, which may be 
responsible for some of the brightness in between the lens galaxies 
in Fig.Hl 

We now explore more extensively this type of binary system. 
We note that the galaxies don't need to have any ellipticity or rel- 
ative mis-alignment in order to have an asymmetric potential. We 
define a simple binary lens model consisting of two identical non- 
singular isothermal spheres of cr = 250 km s^^ and core radii of 
100 pc, Zd = 0.6 and Zs = 1.2. 



6.1.2 Critical curves, caustics and image configurations 

In a recent paper, IShin & Evans! i2008ah studied strong lensing 
by binary galaxies, noting the appearance of several catastrophes. 
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Figure 17. Top: mass distribution of our model of SL2SJ 1405+5502, and 
surface brightness of the predicted images. Bottom: critical curves and caus- 
tics 



Even with isothermal spheres instead of the ellipsoids used in mod- 
elling SL2S J 1405+5502 above, the binary system can produce 
an i nner astroid caustic with two "deltoid" (three-cusped) caus- 
tics dSuvu & Blandfordll2006ll . The latter characterize an elliptic 
umbilic metamorphosis, and can be seen in Fig. [19] as the small 
source-plane features in the right-hand panel, lying approximately 
along the Si -axis. The general caustic structure is the result of the 
mer ging of the two astroid caustics belonging to the two galaxies. 

Ishin&Evanl ( l2008j) studied the different metamorphoses 
undergone as the two mass components are brought closer together. 
In order to give a further illustration of the caustic structure, we 
show in Fig. [TS] the evolution with source redshift of the caustics 
behind our binary model system for a given separation AS = 1.5". 
Plots like this emphasise the three-dimensional nature of the caus- 
tic structure, slices of which are referred to as "caustics" in the 
gravitational lensing literature. This is equivalent to the evolution 
with separation, since the important parameters are the ratio of the 
Einstein radii to the separation and the Einstein radius in- 

creases with the redshift. The condition to obtain an elliptic umbilic 
catastrophe with two isothermal spheres is given bv lShin & EvansI 
gOOSa) . 



M - 



|2 

E,2- 



(8) 



At low source redshift, we observe a six-cusp astroid caustic 
(the aftermath of a beak-to-beak calamity at the merging of the two 
four-cusp caustics). As the source redshift increases, we observe 
another two beak-to-beak calamity (the system is symmetric) lead- 
ing to two deltoid caustics, with the beak-to-beak calamity again 
marking a transition involving the gain of two cusps. These deltoid 
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Figure 19. An "elliptic umbilic image configuration." 

caustics then shrink two points, the two elliptic umbilic catastro- 
phes visible at around Zs = 1.6. Beyond this source redshift we 
then observe two different deltoid caustics (with the elliptic um- 
bilic catastrophe marking the transition between one deltoid and 
the next). For clarity, we have not represented the outer caustics 
until the deltoid caustics meet these curves to form a single outer 
caustic and a lips with two butterflies. In Fig.[T9]we present an "el- 
liptic umbilic image configuration," where four images merge in a 
Y-shaped feature - yet to be observed in a lens system. 

The abundance of such image configurations can be estimated 
very roughly from their fractional cross-section. We estimate ge- 
ometrically the fractional cross-section of the deltoid caustics in 
our model system, compared to the multiple imaging area, as being 
around 10""^. 

6.1.3 Magnification 

In order to study the magnification provided by such binary lenses, 
we again plot the fractional cross-section area with magnification 
above some threshold, against source redshift. The reference cross- 
section was taken to be that where the magnification is ^ 3, which 
corresponds roughly to the multiple imaging area. We make this 
plot for our three fiducial sources. We first note that the cross- 
section for a magnification of 10 is quite big (~0.2): here is a high 
probability of attaining a high total magnification with such lenses. 
This probability persists even with large source sizes, although the 
cross-section for even higher magnification does decrease quite 
quickly. 

For small sources, = 0.5 kpc, the pair of two cusps meeting 
in a beak-to-beak metamorphosis dominates the magnification ~ 
10 cross-section, with the maximum cross-sectional area occurring 
when the two cups are close but before the beak-to-beak calamity 
itself. Re-plotting the same small source cross-section with higher 
magnification thresholds (25, 30, 34, 40, Figs. 1201 and 12 it. we see 
that the magnification ~ 30 — 40 cross-section is dominated by 
the deltoids and their neighbouring cusps. To reach even higher 
magnifications requires a very high redshift (and hence small in 
solid angle) source placed at the center of the main astroid. 

6.1.4 Two merging galaxies 

As an extreme case of the binary system illustrated above, we 
briefly consider two merging galaxies. The model used is simi- 
lar except that we use elliptical profiles with ej, = 0.2 instead of 
spheres. 

In Fig. (22] we plot the convergence and an enlarged view of 
the caustics for a mass component separation of = 0.4 arc- 
seconds and a relative orientation (j) = 1-55. When we vary the 




redshift 

Figure 21. As Fig.|20] but showing higher magnification thi'esholds. 




Figure 22. Two elliptical galaxies with a fixed separation = 0.4 arcsec- 
onds and <j> = 1-55 radians. Lens convergence map left, and a zoomed-in 
view of the astroid caustic in the source plane right. 

orientation, the fold breaks leading to a swallowtail, then a butter- 
fly and then a region of nine-images multiplicity (as presented in 
the Fig.l22t. similar to that seen in the disk and bulge model Sec- 
tion l5.ll In Fig.[23] we then plot the convergence, the critical curves 
and the caustics for zero separation and a relative orientation of 
<j) = 1.2, as might be seen in a line-of-sight merger or projection. 
This model has boxy isodensity contours (as also seen in the crit- 
ical curves), while the caustics include two small butterflies. This 
model is therefore a a physical way to obtain the box i ness needed to 
genera te such catastrophes (e.g. lEvans & Witll200ll:IShin & EvansI 
l2008al) . 

6.2 A four-component lens system: SL2SJ0859— 0345 

6.2.1 Model 

SL2SJ0859— 0345 is an example of a lens with still more complex 
critical curves and caustics; we represent it qualitatively with the 
following four-component model: 



type 






rc 


Oi 


92 


<^ 


MS 


250 


0.0 


0.1 


-1.3172 


0.87 





NIE 


330 


0.1 


0.1 


1.9458 


0.12 


1.2 


NIE 


200 


0.1 


0.1 


-1.2872 


-2.17 


-0.3 


MS 


120 


0.0 


0.1 


0.6885 


1.17 






where once again a is in km s^^, rc in kpc, 6 in arcseconds and cj) 
in radians. The redshifts chosen were Zd = 0.6 and Za ~ 2. The 
velocity dispersions and the ellipticities were chosen to reflect the 
observed surface brightness. In Fig.|24]we show our best predicted 
images, for comparison with the real system in Fig.|9] 
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Figure 18. Evolution of the caustics with source redshift behind a fixed-separation binary lens system. 
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Figure 20. Fractional cross-section for magnification higher than a given threshold (compared to area with ^ 3) for the binary lens model. We consider 3 
fiducial source sizes, from left to right: rc = 0.5, 1.5, 2.5 kpc. 




eral beak-to-beak calamities result in the extended and folded-over 
central astroid caustic, while the elliptic umbilic catastrophe oc- 
curs at around Zs — 2.3. As a caveat, this four-lens model is in one 
sense maximally complex: any smooth group-scale halo of matter 
would act to smooth the potential and lead to somewhat simpler 
critical curves and caustics. However, it serves to illustrate the pos- 
sibilities, and may not be far from a realistic representation of such 
compact, forming groups. 




Figure 24. Convergence and predicted images for a qualitative model of 
SL2SJ0859— 0345, showing the 4 lens components and complex oval sys- 
tem of arcs. 



6.2.2 Critical curves, caustics and image configurations 

In Fig. |25] we present an array of different image configurations 
produced by varying the source position, illustrating what might 
have been. Several different multiplicities are possible, including an 
elliptic umbilic image configuration. The caustic structure is com- 
posed of two outer ovoids, a folded astroid caustic, a lips, and a 
deltoid caustic. The highest image multiplicity obtained is nine (in 
practice eight images, and referred to as an octuplet). We also plot 
the evolution of these caustics with source redshift in Fig.|26l sev- 



6.2.3 Magnification 

We choose to not discuss extensively the magnification of this sys- 
tem, having already examined the available magnifications in the 
simpler instances of the various caustic structures in previous sec- 
tions. Nevertheless, we would still like to know if the deltoid caus- 
tic (when we are close to an elliptic umbilic catastrophe) can lead 
to a higher magnification than other regions of the source plane. In 
order to do that, we compute the total magnification for two partic- 
ular positions of the source: in the central astroid caustic (leading to 
eight observable images in a broken Einstein ring, bottom left panel 
of Fig.|25t, and in the deltoid caustic (top right panel of Fig. [25}. 

We again consider our three fiducial exponential sources (re = 
{0.5, 1.5, 2.5} kpc. The central position leads to magnifications of 
H ~ {55, 39, 27} respectively, and for the deltoid position we ob- 
tain a total magnification of /i ~ {49, 33, 24}. The deltoid caustic 
again does not give the maximum total magnification - it is a re- 
gion of image multiplicity seven, while the central astroid caustic 
of multiplicity nine. However, the deltoid does concentrate, as a 
cusp does, the magnification into a single merging image object. 
In the next section we will see an example of how such focusing 
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Figure 23. Two elliptical galaxies with a relative orientation of 1.2 radians, separation = arcsec. 



by a higher order catastrophe leads to a highly locally-informative 
image configuration. 



6.3 Abell 1703 

Fig.|27]shows the central part of the c luster Abell 1703, and it s criti- 
cal curves and caustics as modeled bylLimousin et al. ' ('2008b'). This 
elliptical cluster, modelled most recently by Oguri et al. ( 2009) and 
?, contains an interesting wide- separation (~ 8 arcsec) 4-image 
system (labelled 1.1-1.4 in the left-hand panel of Fig. |27|(. This 
"central ring" quad has a counter-image on the opposite side of the 
cluster (labelled 1.5). Two red galaxies lie inside the ring; they are 
marked as two red dots in the central panel of Fig.|27| 

Explaining the central ring as one image of a double system 
being sp lit by the two red galaxies is difficult; in the lLimousin et al.] 
( l2008bl) mass model they are represented as isothermal spheres 
with velocity dispersion ~ 100 km s~^ each. The wide separa- 
tion of the ring is a feature of the cluster-scale mass distribution, 
as we now show. We made a rough modeling of the cluster using 
a NIE profile with velocity dispersion of ~ 1200 km s~^, a core 
radius of ~ 30 kpc, and an ellipticity of the potential of ej, = 0.2. 
This single halo model reproduces the central ring image config- 
uration very well. We show the predicted images, critical curves 
and caustics in the upper panels of Fig. |28] This image configura- 
tion is associated with a hyperbolic umbilic catastrophe; one can 
see in the right-most panel of Fig.[27]and more clearly in he corre- 
sponding panel of Fig.|28]that the source lies very close to the point 
where the ovoid and astroid caustics almost touch. Were the source 
actually at the catastrophe point, the four images of the central ring 
would be merging; as it is, we see the characteristic pattern of a 
hyperbolic umbilic image configuration. As far as we know this is 
the first to have ever been observed. 

When we add the two a — 100km s~^ SIE galaxies inside the 
central ring (lower panels of Fig. [28}, th e central ring system is p er- 
turbed a little. This perturbation allowed Limo usin et alj j2008bh to 
measure them despite their small mass, illustrating the notion that 
sources near higher order catastrophes provide opportunities to map 
the local lens mass distribution in some detail. Indeed, this image 
system provides quite a strong constraint on the mass profile of the 
cluster dLimousin et al.|[2008bl) . 



6.3.1 An even simpler model 

To try and estimate the abundances of such hyperbolic umbilic im- 
age configurations, we consider a slightly more representative (but 
still very simple) model of a cluster: cluster and its brightest clus- 
ter galaxy (BCG), both modelled by a NIE pro file. The parameters 
used are taken after considering the sample of ISmith et alj j2005l) 
who modelled a sample of clusters using NIE profiles. The model 
is as follows: 



type 




(7 


/ 




NIE- 


Cluster 


900 


0.2 


80 


NIE- 


BCG 


300 


0.2 


1 



where a is in km s and Tc is in kpc. 

6.3.2 Critical curves, caustics and image configurations 

Using a simple NIE profile, we can obtain an hyperbolic umbilic 
catastrophe, provided that the central convergence is sufficiently 
shallow ( i.e. a significant core r adius is required), and sufficiently 
elliptical. iKormann et al.l ( Il994h give a criterion for a NIE lens to 
be capable of generating a hyperbolic umbilic catastrophe: 



where 6c is the core radius in units of the Einstein radius, and the 
axis ratio q refers to the m ass distribution and not the lens potential. 
Kassi ola & Kovneil ( Il993h also discuss elliptical mass distributions 
and elliptical potentials, and describe the different metamorphoses 
when varying the core radius. 

We show the evolution of this model's caustics with source 
redshift in Fig. [29] At low redshift, the caustic structure consists of 
two lips (with relative orientation of 90 degrees); as Zs increases, 
the inner lips becomes bigger and meets the outer lips in a hyper- 
bolic umbilic catastrophe. The hyperbolic umbilic metamorphosis 
conserves the number of cusps: at higher source redshift, we see the 
familiar (from elliptical galaxies) outer ovoid caustic with an inner 
astroid caustic. 

In the light of these considerations, we don't expect hyperbolic 
umbilic catastrophes to be generic in elliptical galaxy lenses, as 
the central convergence is not sufficiently shallow. However, this 
catastrophe point should be present behind many clusters, as the 
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Figure 25. Alternative histories for tlie SL2SJ0859— 0345 system: the image configurations (left) and multiplicities resulting from various source positions 
(right). Left column, from top to bottom: 1 . Double image, of a source in the first outer caustic. 2. Quad: second outer caustic. 3. Sextuplet: outer astroid caustic. 
4. Quad: lips caustic. 5. Octuplet: inner-astroid caustic. Right column, from top to bottom: 1. Sextuplet, including 4 merging images of an elliptic umbiUc 
configuration: source in a deltoid caustic. 2. Octuplet: source at the astroid fold crossing point. 3, 4, 5. Octuplets: sources on each of the three inner-astroid 
cusps. 
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Figure 26. Evolution of the caustics witli source redshift for our qualitative model of SL2SJ0859— 0345. The ovoid outer curves are sometimes only sparsely 
sampled by the adaptive mesh code used to compute the magnification maps (see appendix). 
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Figure 27. Abell 1703 and the lens model modeled of lLimousin et alj fcOOSbh Left: the core of the cluster obsei-yed wit h HST/ACS jStotlll2007l) . The image 
is approximately 40 arcsec wide. Critical curves (centre) and caustics (right) in the lens model o f iLimousin etal] j2008bl) . The two galaxies inside the central 
ring are represented by two filled red circles. 



combination of a BCG and a cluster with appropriate core radii 
and ellipticities can lead to the central convergence and ellipticity 
required. Of course the hyperbolic umbilic image configuration is 
more likely to happen than the catastrophe itself. Indeed, as Fig.|29l 
shows, the hyperbolic umbilic image configuration can occur over 
a wide range of source redshift. 



6.3.3 Magnification 

Fig. [30] shows the fractional cross-section area of our NIE cluster 
model (at Zd = 0.28) as a function of magnification threshold. 
Here, the area of the source plane providing total magnification 
greater than jj, is compared to the area providing total magnification 
greater than fi 15, which is approximately the area of the astroid 
caustic. We consider again our fiducial exponential sources, and 
observe that the system produces some very high magnifications. 
Analyzing the different convolved maps of source plane magnifi- 
cation, we find that it is indeed the regions close to the hyperbolic 
umbilic catastrophes that provide the highest magnifications. For 
rc = 0.5 kpc, the cross-section for /x > 100 is entirely due to 
the two regions close to the catastrophe. We obtain the same re- 
sult for the two other sources (respectively /i ~ 75 and /i ~ 25 
for Tc = 1.5 and 2.5 kpc): it is likely that A 1703 is magnifying 
its source by something like these factors. Fig. [30] shows that the 
region with a characteristic magnification of a hyperbolic umbilic 
catastrophe is around 1-2% of the multiple-imaging area. 




Figure 30. Fractional cross-section (relative to the area with magnification 
/i ^ 15) as a function of the threshold magnification for our simple NIE 
cluster model. 



7 THE ABUNDANCE OF HIGHER-ORDER 
CATASTROPHES 

In the previous sections, we have illustrated a range of "exotic" 
lenses, systems presenting higher-order catastrophes such as swal- 
lowtails, butterflies, elliptic umbilics and hyperbolic umbilics, and 
their related metamorphoses. In this section we make the first very 
rough approximations of the abundances of lenses showing such 
exotic image configurations. 

Procedure. For each given type of lens, we first estimate the 
basic abundances of the more common, lower magnification config- 
urations, based on the number densities observed so far in SLAGS 
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Figure 28. Simple models for understanding Abell 1703, with predicted images (left), critical curves (centre) and caustics (right). The effect of the two galaxies 
inside the central ring is shown in the lower panels; they are absent from the model in the upper panels. 
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Figure 29. Evolution of a simple elliptical cluster lens caustics with source redshift. The closeness of the cusps to the outer ovoid caustic persists over a wide 
range of redshifts. 



and SL2S. We then estimate, for a given source redshift, the frac- 
tional cross-section area for production of the desired image con- 
figuration. To do this we consider both a) the information provided 
by the magnification maps and their associated fractional cross- 
sectional areas, and b) our rough estimates of the "area of influ- 
ence" of the relevant caustic feature. Finally, we estimate the typ- 
ical window in source redshift for which the image configuration 
is obtainable, and hence the fraction of sources lying at a suitable 
position along the line of sight. The product of these two fractions 
with the basic abundance gives us an order-of-magnitude estimate 



of the number density of exotic lenses detectable in a future imag- 
ing survey (which we take to cover ~ 10* square degrees). 

Galaxy scales. On galaxy scales, we expect some 10 strong 
gravitational lenses observable per square degree, provided the im- 
age resolution is hig h enough: this is approximat ely what is found 
in the COSMOS dFaure et alj boOSh and EGS jMoustakas elal] 
2007) surveys. The m ajority of these len s galaxies are massive 
ellipticals, as expected jTumer et al]|l984h . In the larger SLACS 
sample of similar galaxy-galaxy lenses, the fraction of systems 
showing strong evidence of a disk component is approximately 
10%, while only one or two (a few percent) have massive satel- 
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lites within tlie Einstein radius jBolton et al.l2008h . Making the link 
between galaxy-scale lenses expected in imaging surveys and the 
lenses found by the SLACS survey should be done with care, but 
it is not unreasonable to start from this extrapolation - especially 
in the light of the SLACS studies demonstrating the normality of 
the lens galaxies relative to the parent massive galaxy population 
jTreu et aljioo^ : ?). This gives us basic abundances of disky lenses 
and binary/merging lenses of 1 and ~ 0.1 per square degree re- 
spectively. This is what we might expect a space-based survey to 
recover; the number visible from the ground will be much reduced 
by the image resolution. 

Therefore, our estimate of the abundance of butterflies is 
1 X (10/180) X 0.001 ~ 5 X 10"^ per squai-e degree, where 10/180 
is the angular range where we can expect butterflies, 0.001 is our es- 
timation of the cross-section of the butterflies and nine-image mul- 
tiplicity regions from Section lSTI Note that we have optimistically 
considered a uniform distribution of disk-bulge misalignments: this 
abundance might be taken as an upper limit for these types of lens 
mass distributions. We expect to have roughly 10 times more swal- 
lowtails, since the range of misalignment angles leading to swal- 
lowtails is ~ 10 times greater than that for butterflies. We therefore 
obtain an upper limit of ~ 1 observable disk/bulge butterfly on the 
sky, and ~10 disk/bulge swallowtails. 

Group scales. On group scales, we might expect the ground- 
based surveys to be as complete as those from space. Moreover, the 
depth of the CFHT Legacy Survey images is comparable to that ex- 
pected from future surveys such as that planned with LSST (which 
will cover 2 orders of magnitude more sky area). Extrapolating the 
number of group- sca le gravitational lenses from the SL2S survey 
jCabanac et ai]|2007h is therefore a reasonable starting point. The 
observed abundance of group-scale lenses is ~ 0.1 per square de- 
gree, and perhaps a third show multiple mass components within 
the Einstein radius. To order of magnitude we therefore estimate 
the basic abundances of compact group lenses as ~ 0.01 per square 
degree. 

The two systems studied here show a variety of complex caus- 
tic structure. In the case of the binary systems in Section [6Tl we es- 
timated geometrically the fractional cross-section of a deltoid caus- 
tic to be ~ lO^"* over a wide redshift range. We therefore obtain an 
order of magnitude estimate of 0.01 x 10^'^ ~10^^ per square de- 
gree for the abundance of elliptic umbilic image configurations, or 
approximately one observable Y-shaped image feature in our future 
survey. This cross-section area is unlikely to be very different (de- 
pending on the lens) to that for the production of swallowtails and 
butterflies by group lenses. A Monte Carlo simulation would per- 
haps be the best way to estimate these abundances more accurately. 

Clusters. For clusters, as shown in section lF.3.21 the properties 
needed to obtain an hyperbolic umbilic catastrophe are well under- 
stood: the problem lies in estimating the ranges of ellipticities, con- 
vergences and core radii of the cluster and BCG mass distributions 
that combine to provide a sufficiently shallow inner density slope 
for the system to generate a hyperbolic umbilic. As a first attempt, 
we consider the distribution of cluster lens potential ellipticities es- 
timated from numerical simulations: this is appro ximate ly a Gaus- 
sian with mean 0.125 and width 0.05 ( Marshaii|[200l following 
Ijing & Sutoll2002h . In general the BCG increases the central mass 
slope (decreasing the core radius of an effective NIE profiles), such 
that we need high ellipticity in order to be close to an hyperbolic 
umbilic catastrophe. We estimate that only the 1-10% highest el- 
lipticity clusters will meet this criterion. Our qualitative analysis 
in Section 16.31 suggests that once the ellipticity and density pro- 
file condition is met then the hyperbolic umbilic metamorphosis 



is slow with source redshift. The fractional cross-section area of 
this "hyperbolic umbilic region" is approximately 2%. As the num- 
ber of cluster lenses is expected to be ~ 0.1 per square degree, 
we estimate the abundance to be ~ 10^* per square degree, again 
predicting ~ 1 hyperbolic umbilic catastrophe image configuration 
per sky survey. A question comes to mind: is Abell 1703 the only 
hyperbolic umbilic gravitational lens we will ever see? 

Discussion. Throughout this work we have focused on the 
most abundant sources, the faint blue galaxies, and tried to ac- 
count for their finite size when estimating magnifications and 
now abundances. To first order this should take care of any 
magnification bias. However, for the less common point-like 
sources (such as quasars and AGN) it is possible that magni- 
fication bias may play a much more important role. Investiga- 
tion of this effect would be a good topic for further work. The 
sensitivity of galaxy-scale gravitational tensing of p oint sources 
to smril-scale CDM substructure (e.^. Mao & Sch neider 199^; 
iDalal & Kochanekll2002l : lBradac et alj2004ah is such that more in- 
vestigation of the exotically-lensed quasar abundance is certainly 
warranted. There are similar, but perhaps more constraining, rela- 
tions between the signed fluxes of images of sources lying close to 
higher order cata strophes ( Aazami & Fetters 2009) as there are for 
folds and cusps HKeeton et all2003l.l2005h . 

Another possibility that might be included in a more com- 
plete analysis of exotic lens abundance is that of generating higher 
order catastrophes with multiple lens planes, a topic studied by 
iKochanek & Apostolakis 1988). If a substantial fraction of the ob- 
served group-scale lenses are in fact due to superpositions of mas- 
sive galaxies rather than compact groups this could prove a more 
efficient mechanism for the production of higher order catastro- 
phes. 



8 CONCLUSIONS 

Inspired by various samples of complex lenses observed, and with 
future all-sky imaging surveys in mind, we have compiled an atlas 
of realistic physical gravitational lens models capable of producing 
exotic image configurations associated with higher-order catastro- 
phes in the lens map. For each type of lens considered, we have 
investigated the caustic structure, magnification map, and exam- 
ple image configurations, and estimated approximate relative cross- 
sections and abundances. We draw the following conclusions: 

• Misalignment of the disk and bulge components of elliptical 
galaxies gives rise to swallowtail (for a wide range of misalignment 
angles) and butterfly catastrophes (if the misalignment is close to 
90 degrees). The image configuration produced by the butterfly 
caustic would be observed as a broken (8-image) Einstein ring. 

• The central nine-imaging region has cross-section '--^ 10^'^ rel- 
ative to the total cross-section for multiple imaging; combining this 
with rough estimates of the abundance of such disky lens galax- 
ies in the SLACS survey, we estimate an approximate abundance 
of < 1 such butterfly lens per all-sky survey. We estimate that the 
swallowtail lenses could be ~ 10 times more numerous. 

• Binary and merging galaxies produce elliptic umbilic catastro- 
phes when the separations of the mass components are comparable 
to their Einstein radii. The configuration formed when the source 
lies within the deltoid caustic is a Y-shaped pattern of 4 merging 
images between, but offset from, the two lenses. 

• In this case the deltoid caustic does not provide the maximum 
total magnification available to the lens; the relevant cross-section 
must be estimated from the area within the deltoid. We find the 
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relative cross-section to also be ^ 10~ , leading to an approximate 
abundance of ~ 1 binary galaxy elliptic umbiUc lens per all-sky 
survey. 

• More complex group-scale lenses, of the kind being discov- 
ered by the SL2S survey, offer a wide range of critical points in 
their caustics, and so are a promising source of exotic lenses. For 
example, our model of SL2SJ0859— 0345 shows an elUptic umbiUc 
catastrophe at around z = 2.4. At present, our understanding of the 
distribution of group-scale lens parameters, and the extreme vari- 
ety of caustic structure prevents us making accurate estimates of 
the exotic lens abundance for these systems. 

• As noted by previous authors, elliptical clusters with appropri- 
ate inner density profile slopes produce hyperbolic umbilic catas- 
trophes. We find that just such a simple model is capable of repro- 
ducing the central ring image configuration of Abell 1703, demon- 
strating the source to lie close to a hyperbolic umbilic point. The 
total magnification of this source could be ~ 100, depending on 
the source size. 

• Despite the rather general properties of clusters needed to pro- 
duce such image configurations, we estimate that the all-sky abun- 
dance of such hyperbolic umbilic cluster lenses may still only be 
~ 1. 

In some cases, proximity to a catastrophe point does not guar- 
antee maximal magnification. Further studies of the abundances of 
exotic lenses should explore more advanced diagnostics of catas- 
trophic behavior. While high magnification is certainly a desirable 
feature of gravitational lenses for some applications (e.g. cosmic 
telescope astronomy), high local image multiplicity is perhaps the 
more relevant property for studies of small scale structure in the 
lens potential, using the information on the gradient and curvature 
of the mass distribution that is present. The constraining power of 
the hyperbolic umbilic image configuration in Abell 1703 is a good 
example of this on cluster scales. The exploitation of the plausible 
lenses we have predicted here in this way would be an interesting 
topic for further research. 
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APPENDIX A: LENS AND SOURCE MODELS 

In this appendix, we describe in a litfle more detail the different 
models used throughout the paper: the isothermal profiles to model 
the lens potentials and the Sersic profiles for the background source 
galaxies' light distributions (or for the stellar mass component of a 
lens galaxy). 

For all image simulations in this paper we uses the publically- 
available GLAMROC code (Gravitational Lens Adaptive Mesh Ray- 
tracing of Catastrophes) written by E. A. Baltz|j Here we give a 
very brief introduction to this code and its capabilities. 

Each GLAMROC lens model is made of an arbitrary number 
of lensing objects on an arbitrary number of lens planes. The in- 
dividual lensing object all have analytic lens potentials, such that 
all their derivatives are also analytic. In this way the deflection an- 
gles, magnification matrices, and combinations of higher deriva- 
tives can be calculated as sums of terms coming from each lens in 
turn. Higher derivatives are required to identify catastrophes of the 
lens map, which include all the catastrophes define above. An adap- 
tive mesh is used to improve the resolution either near the critical 
curves of the lens system or in regions of high surface brightness. 
Lens types currently implemented include point masses, isother- 
mal spheres with and without core and truncation radii, NEW pro- 
files with and without truncation, and Sersic profiles. Eor each type, 
elliptical isopotentials can be used and either boxiness, diskiness 
and skew added. Sources are modelled with (superpositions of) 
elliptically-symmetric Sersic profiles. 



Al Isothermal lenses 

Throughout the paper, we used the non-singular isothermal el- 
lipsoid (NIE) model for our lens components. The NIE (three- 



dimensional) mass density diverges as p cc r 
jected) mass density profile is 

E(r) 



2G(r2 +r2)i/2^ 



'; the surface (pro- 



(Al) 



where r is now the projected radius, rc is the core radius, and a 
is the velocity dispersion of the lens. The core radius is that at 
which the rising density profile turns over into a uniform distri- 
bution in the central region. This type of profile has been shown to 
be a very good approximation of the lens potential on both galaxy 
scales (e.g. lKoopmans et al.l200 (/). provided the core radius is very 
small (< 0.1 kpc), and als o on cluster scale s, with much larger 
core radii (~ 50.0 kpc, e.g. ISmith et aT I l2005h . In practice a small 
core allows more convenient plotting of caustics and critical curves. 
However, we are careful not to allow cores larger than typically 
permitted by the data: the non-singular core prevents the infinite 
de-magnification of the central image, which would be observable. 
The NIE model isopotential lines are ellipses with constant 
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where x and y denote the Cartesian coordinates in the lens plane. 
The potential is, therefore, a function of r only, and the ellipticity e'p 
is the ellipticity of this potential. The ellipticity must be carefully 
chosen in order to keep the potential physically meaningful. Indeed, 
for e'„ > 0.2 the isodensity co ntours become dumbbell-shaped (see 
e.g. iKassiola & Kovneilll993l) . ? provide a simple practical solu- 
tion to this problem: add several (suitably chosen and weighted) 
elliptical potentials at the same location. Using their algorithm, we 



are able to model nearly elliptical isodensity contours with ellip- 
ticities as large as e' = 0.8. This allows us to model, for ex- 
ample, the combination of a bulge and a disk in a galaxy. (Note 
th at the NIE model desc ribed here is different from that defined 
bv lKormann et alj ( ll994l) . where the ellipticity pertains to the mass 
distribution, not to the lens potential.) 

We follow the GLAMROC notation and use as our main ellip- 
ticity definition e — f^npp'' where a and b are the major and minor 
axis lengths of the ellipse in question. This is different from the def- 
inition often used in weak lensing, e = . The relation between 



these two definitions of the ellipticity is e 



_2e 

1 + 



A2 Sersic sources 

The Sersic power law profile is one of the most frequently used 
in the stud y of galaxy mor phology. It has the following functional 
form (e.g. iPeng et al.[|2002[) . 
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where Ic is the surface brightness at the effective radius rc. The 
parameter is known as the effective, or half-light, radius, defined 
such that half of the total flux lies within rc. The parameter n is the 
Sersic index: with n = 1/2 the profile is a Gaussian, n = 1 gives 
an exponential profile typical of elliptical galaxies, and n = 4 is 
the de Vaucouleurs profile found to well-represent galaxy bulges. 
We compared these three standard profiles when exploring the ob- 
servability of exotic lenses in the text. 

The half-light radius is an important parameter when consider- 
ing the observable magnification of extended sources. For the sizes 
of the expected faint blue s ource galaxies, we ad opted the size- 
redshift relation measured bv lEerguson et al.l ( |2004|) . For reference, 
this gave effective radii of 4.0 kpc (0.5") at Zs = 1, decreasing to 
2.9 kpc (0.35") at = 2 and 1.7 kpc (0.21") at = 4. We used 
sources having equal half-light radii, instead of the same total flux 
(the relative intensity is, as we explain in the paper, not important 
as the magnification calculation based on Equation|7]is independent 
of the total flux of the source). 



^ http : / /kipac . st an ford. edu/ col lab/ research/ lensing/ glamroc/ 



